Universal low-temperature tricritical point in metallic ferromagnets and ferrimagnets 
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An earlier theory of the quantum phase transition in metallic ferromagnets is revisited and gener- 
ahzed in three ways. It is shown that the mechanism that leads to a fluctuation-induced first-order 
transition in metallic ferromagnets with a low Curie temperature is valid, (1) irrespective of whether 
the magnetic moments are supplied by the conduction electrons or by electrons in another band, (2) 
for ferromagnets in the XY and Ising universality classes as well as for Heisenberg ferromagnets, and 
(3) for ferrimagnets as well as for ferromagnets. This vastly expands the class of materials for which 
a first-order transition at low temperatures is expected, and it explains why strongly anisotropic 
ferromagnets, such as UGe2, display a first-order transition as well as Heisenberg magnets. 
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I. INTRODUCTION, AND RESULTS 

Quantum phase transitions are a subject of great 
interest .Gill In contrast to classical or thermal phase tran- 
sitions, which occur at a nonzero temperature Tc > and 
are driven by thermal fluctuations, quantum phase tran- 
sitions occur at zero temperature, T = 0, as a function 
of some non-thermal control parameter and are driven 
by quantum fluctuations. In this paper we will focus 
on quantum phase transitions in metallic systems. For 
reasons discussed below, these transitions are especially 
interesting. 

A prototypical quantum phase transition is the one 
from a paramagnetic metal to a ferromagnetic metal. 
Indeed, the earliest theory of a quantum phase transi- 
tion was the Stoner theory of ferromagnetism.^ Stoner 
assumed that the conduction electrons are responsible 
for the ferromagnetism, and developed a mean-field the- 
ory that describes both the classical and the quantum 
ferromagnetic transition. In an important paper. Hertz 
later derived a Landau-Ginzburg-Wilson (LGW) func- 
tional for this transition by considering a simple model 
of itinerant electrons that interact only via a contact po- 
tential in the particle- hole spin-triplet channel.^ Hertz 
analyzed this (dynamical) LGW functional by means of 
renormalization-group (RG) methods. He concluded that 
the critical behavior in the physical dimensions d — 2 and 
d = 3 is mean-field-like. That is, as far as the static crit- 
ical exponents of the transition at T = are concerned, 
he concluded that Stoner theory is exact in d = 2 and 
d = 3. 

In the mid 1990s it was realized that the above con- 
clusion is not correct. The problem is that in metals 
at T = there are gapless particle-hole excitations that 
couple to the magnetic order-parameter fluctuations and 



influence the quantum critical behavior for all dimensions 
d < 3. In Hertz's theory this coupling is taken into ac- 
count only in an approximation that does not suffice for 
yielding the leading critical behavior. Technically, Hertz 
theory treats the fermionic soft modes in a tree approx- 
imation, whereas describing their influence on the criti- 
cal behavior requires taking into account fermionic loops. 
Physically, a correct description of any phase transition 
must treat the order parameter fluctuations and all soft 
modes that couple to them on equal footing. 

A theory that takes into account these effects was de- 
veloped by the present authors and T. Vojta. In Ref. |4| 
it was shown that the quantum phase transition from a 
metallic paramagnet to an itinerant ferromagnet in the 
absence of quenched disorder in d = 2 and d = 3 is 
generically discontinuous, or of first order, in contrast to 
the second-order transition with mean-field critical be- 
havior predicted by Hertz theory.lSl The mechanism be- 
hind this phenomenon is analogous to what is known 
as a fluctuation-induced first-order transition in super- 
conductors and liquid crystals.^ There, soft fluctuations 
of the electromagnetic vector potential (in superconduc- 
tors) or the nematic order parameter (in liquid crystals) 
couple to the order parameter and effectively change the 
sign of the cubic term in the equation of state, leading to 
a first-order transition. In the quantum magnetic case, 
the role of the additional soft modes is played by the 
fermionic particle-hole excitations mentioned above that 
are massless at T = 0. Since these modes acquire a mass 
at T > 0, the tendency towards a first-order transition 
diminishes with increasing temperature. This leads to a 
tricritical point at a temperature T^c > that separates 
a line of continuous transitions at T > Ttc from a line 
of first-order transitions at T < Ttc- In a later paper 
with RoUbiihler, the effects of a magnetic field H were 
investigated.— It was found that in the space spanned by 
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FIG. 1: Generic phase diagram of a metallic magnet in the 
space spanned by temperature {T), magnetic field (H), and 
the control parameter (i). Shown are the long-range or- 
dered magnetic (LRO) and paramagnetic (PM) phases, lines 
of second-order transitions, surfaces of first-order transitions 
( "tricritical wings"), the tricritical point (TCP), and the two 
quantum critical points (QCP). The long-range order can be 
of ferromagnetic or ferrimagnetic type, and the electrons caus- 
ing the long-range order can be in the same band as the con- 
duction electrons, or in a different band. See the text for 
further explanation. 



[TJalso applies to generic metallic ferrimagnets. Ferrimag- 
nets are materials that spontaneously develop both a ho- 
mogeneous and a staggered magnetization at the same 
critical value of either the temperature (for a classical 
transition) or a non-thermal control parameter (for a 
quantum transition). Physically, this can happen when 
magnetic moments of unequal magnitude on a bipartite 
lattice align in opposite directions.^ 

The unifying principle behind these generalizations is 
the realization that coupling a homogeneous magneti- 
zation to conduction electrons will produce the same 
results irrespective of the microscopic origin of the 
magetization.HS As a result, the phase diagram depicted 
schematically in Fig. [T] is valid for generic metallic ferro- 
magnets in addition to itinerant ones, for ferromagnets of 
XY or Ising type in addition to Heisenbcrg magnets, and 
for ferrimagnets as well as for ferromagnets. In all cases 
we also consider the effects of nonmagnetic quenched dis- 
order. In Ref. ,4| it was shown that this type of disor- 
der leads to an interesting phase diagram with a num- 
ber of multi-critical points, and that sufficiently strong 
quenched disorder causes the first-order paramagnetic- 
to-ferromagnetic transition in metals to become second 
order. We will see that the same result holds for metallic 
ferrimagnets. Experimentally, the effects of disorder on 
either one of these transitions have not yet been studied 
systematically. 



T, H, and the control parameter, tricritical wings, or 
surfaces of first-order transitions, emanate from the tri- 
critical point and terminate in a pair of quantum critical 
points in the T = plane. The wing boundaries at T > 
are given by lines of critical points that are reminiscent 
of a conventional liquid-gas critical point and connect 
the tricritical point with the quantum critical points at 
T = 0. The resulting generic phase diagram is shown in 
Fig. IT] This general picture is in good agreement with 
experimental results for low-Curie-temperature metallic 
ferro magne ts, including ZrZn2,'^ UGe2,'^ URhGe,^" and 
MnSi.EEl 

In this paper we generalize our previous theory in three 
important ways. First, we show that our previous results, 
which had been derived under the same assumption made 
by Stoner and by Hertz, namely, that the magnetism is 
caused only by itinerant electrons, remain valid in metal- 
lic systems where the magnetism is caused by electrons 
in a different band than the conduction electrons. 

Second, we show that the results are not restricted 
to Heisenberg ferromagnets, contrary to what was im- 
plied in Refs. S] and [131 Rather, they apply equally 
well to metallic XY or Ising magnets, since the mag- 
netic moments couple to conduction electrons whose 
spins have three degrees of freedom. This is an im- 
portant point, since some of the relevant materials are 
strongly anisotropic magnets, including UGe2 (easy axis) 
and URhGe (easy plane). 

Third, we show that the phase diagram shown in Fig. 



II. THEORY 

We now derive the results listed in Sec. HI To this end 
we are interested in a theory that describes the magne- 
tization or order-parameter (OP) field M, the fermionic 
degrees of freedom described by Grassmann- valued fields 
■0 and V'l a-nd the coupling between them. Accordingly, 
the action will have three parts: 

A[M;'iP,'tP] ^Aop[M]+AF[i^,^P]+Ac[M;i^,^p] , 

(2.1a) 
and the partition function is given by 



Z = D[M]D[iP,il;] 



-A[M-il>,t(>] 



(2.1b) 



We are, however, not interested in a complete descrip- 
tion of the fermionic degrees of freedom; rather, we want 
to restrict ourselves to the fermionic soft modes and in- 
tegrate out the massive modes in the simplest approxi- 
mation that respects the symmetries of the problem to 
arrive at an effective Landau-Ginzburg-Wilson (LGW) 
theory in terms of soft modes only. If we denote the soft 
fermionic degrees of freedom collectively by q, and the 
massive ones by P, we formally have 



Z = J D[M,q] 



AM,q] 



(2.2a) 



where 



^LGw[M,(j] = Aop[M]-\n f D[P] 



xe 



-^Fk,P] 



-A,[M:q,P] 



= ^op[M]+^FM+A[M,g] .(2.2b) 

As we will see later, the q are matrices forraed by bilinear 
products of the fermion fields, qnm{x, y) = 'tpnix) ipmiv) 
with {n + 1/2) (m + 1/2) < 0, and the P are given by 
the same products with (n + 1/2) (m + 1/2) > 0. Here 
ipnix) = ip{x,ujn) is the temporal Fourier transform of 
the Grassmann field ip{x), where x = {x,t) comprises 
the real-space position x and the imaginary-time variable 
T in a Matsubara formalism, and w„ = 2T:T{n + 1/2) 
is a fermionic Matsubara frequency, ipnix) is defined 
analogously. 

This separation of soft and massive fermionic modes 
q and P, respectively, integrating out P in a suitable 
approximation, and determining the consequences of the 
coupling between q and M, is the central objective of 
this paper. For the separation we will make use of the 
general theory developed in Ref. [TBI 



A. Order parameter, and coupling to fermions 

We are interested in magnetic order, and hence the 
appropriate order-parameter field is the magnetization 
M{x). We write the magnetization as a part m{x) whose 
average is the homogeneous magnetization, and a part 
n(a;) whose average is a staggered magnetization. 



transition is characterized by both m and n acquiring a 
nonzero expectation value at the same point in parame- 
ter space. In this sense there is only one order parameter 
field for a ferrimagnetic transition; this fact will be impor- 
tan later. For the purposes of the present paper, a crucial 
question is the coupling of the order-parameter fluctua- 
tions to the soft fermionic degrees of freedom. Since the 
soft parts of the latter are soft at zero wave number, 
the leading coupling is to m. The fermions also couple 
to n, but this leads to subleading effects since the stag- 
gered magnetization is soft at a nonzero wave number. 
We will neglect this coupling in what follows. Physically, 
the near-homogeneous magnetization flutuations act as 
a magnetic field proportional to m that couples to the 
electronic spin density 



ns{x) 



a.b 



1paix)(Tabtpbix) 



(2.4a) 



Here cr — (cr^, (T^,cr^) = (ct^, cr^, cr"^) denotes the Pauli 
matrices, and a,b ^ (tii) = (+li^l) are spin indices. 
The coupling takes the form of a Zeeman term 



Ac[M-^p,ijj] = c 



/ dx m{x) ■ ns{x) , (2-4b) 



with c a coupling constant. As we will see, the spin den- 
sity contains both massive and massless modes, so only 
part of Eq. ( |2.4b[ ) contributes to A[M, q] in Eq. ( |2.2b |. 
We will discuss this separation next. 



M{x) = m{x) + n{x) 



N 

yjcos(fcj 



(2.3) 



Here the kj are N wave vectors that characterize the 
staggered magnetic order, and both m{x) and n{x) are 
slowly varying in space and time. In particular, their 
Fourier expansions contain only wave numbers that are 
small compared to the norms of the kj . 

In a paramagnetic state the expectation values of m 
and n are both zero. At a transition to a ferromagnetic 
state the expectation value of m becomes nonzero while 
that of n remains zero; at a transition to an antifer- 
romagnetic state the converse is true. A ferrimagnetic 



B. Fermionic soft modes 

In this subsection we separate the massless fermionic 
modes from the massive ones by means of the technical 
apparatus developed in Ref. [TH Here we will quote only 
as much of this formalism as is necessary for the further 
development, see Ref. [121 for additional details. 

The soft fermion excitations are all two-particle ex- 
citations; the related correlation functions are those of 
bilinear products of fermion fields. The latter commute 
with each other, and with individual fermion fields, and 
hence are isomorphic to classical fields. Denoting these 
classical fields by Q, we define a classical matrix field 



^nm\X-, y } 



I -'lljnt{x)frnt{y) 
-1pnl{x)tlj,nt{y) 

1pnl(x)iimt{y) 
\ -■ipntix)lpmtiy) 



-i'ntix)tpmiiy) 
-^nl{x)lpml{y) 

ilMi{x)iljjni{y) 

-1pntix)tpmi{y) 



~-0„^(a;)V'm;(y) 

-tnl{x)ll^ml{y) 

1p_nl{x)^ml{y) 

-^n^{x)^m\.{y) 



1pnt{x)'tlJrnt{y) \ 

ip_niix)i(;^^{y) 

-1pnl{x)lpmt{y) 
'>Pnt{x)lpmt{y) J 



(2.5) 



r 



Here "=" means "isomorphic to" ; technically, the isomor- phism is implemented by means of a Lagrange multiplier 



field, see below. We also define the Fourier transform of 

Qnm{k,p) = ^ J dxdy e-'^'^+'Py Q,,^{x,y) . (2.6a) 
It is further useful to define 



,{k;q) = Q„^{k + q/2,k~q/2) (2.6b) 



and 



q k 

(2.6c) 
The 4x4 matrix Qnm can be expanded in a spin- 
quaternion basis 

3 
Qnm{x,y) = ^ (Tr«'s,);Q„,„(a;,y) , (2.7) 



where tq = sq = I2 is the unit 2x2 matrix, and 
^1,2,3 = — si,2,3 = — icT^'^'"^. An explicit inspection of 
the 16 matrix elements shows that r = 0, 3 represents 
the particle-hole channel, i.e., products of the form f/;-!/), 
whereas r = 1, 2 represents the particle-particle channel, 
i.e., products of the form ipij} or ^pip. For our purposes 
we will need only the particle-hole degrees of freedom. 

It was shown in Ref . Hn] (see also Ref . [TT]) that a crucial 
criterion for separating the fermionic degrees of freedom 
into soft and massive modes is given by the relative signs 
of the frequency arguments of the matrix elements Qnm- 
Accordingly, we write 

IQnmi^) = r9„™(a;)0(-w„a;„)-|-;P„„(a;)e(a;„w,„) 



(z = 1,2,3) 



(2.8) 



Here Q is the step function, and we use the fact that in 
the spin-triplet channel (i — 1,2,3) the expectation value 
of the Q-matrix vanishes (this is since the fermionic de- 
grees of freedom described by Q do not by themselves 
have long-ranged magnetic order; see the discussion at 
the end of the current subsection) , so that q and P rep- 
resent fluctuations. In what follows we will absorb the 
step functions into the matrix fields q and P, i.e., writ- 
ing qnm implies n > and m < Q and Pnm implies either 
n > and m > or n < and m < Q. The \.q are the 
spin-quaternion elements of a matrix 



g™ 



^{X) ^Y^i'^r 



(^Si)lq„^{x) 



(2.9a) 



It is also useful to define an adjoint matrix 

i,r 

where t^ and sf are the hermitian conjugates of r^ and 
Si, respectively. In addition, the theory contains a field 



i'nm{x) that has the same properties as qnm{x) except 
for different propagators, see below. The origin of ^ is 
the Lagrange multiplier field A that constrains the bilin- 
ear products of fermion fields to the q. In various places 
in the theory q — X = ^ appears, and the A-propagator 
equals minus the g-propagator for noninteracting elec- 
trons, whereas cross-correlations between q and A vanish. 
The net effect of A is therefore to subtract the noninter- 
acting part of the g-propagator wherever the combination 
q — X occurs. 

The q correlation functions are the basic soft modes in 
the theory, see below. However, due to nonlinear cou- 
plings the P couple to the q and thus have a soft compo- 
nent. This effect can be expressed by expanding P in a 
power series in q. To quadratic order in q and to lowest 
order in the fermion interaction one finds 

Pi2(fe) w --2i"^"^(p[%{p,k-p)(p^r^{p)ip^^{k~p) 
3 p 

x[ii3{p)it2{k-p)+iUp)h2{k-p)] . (2.10) 

Here and it what follows we use a simplified notation for 
frequency indices, 1 = ni, etc. We have dropped con- 
tributions to P of higher order in g, and a contribution 
that is linear in the interaction and linear in q, see Ref. 
[TCI neither will be needed for our purposes. We also have 
omitted a term quadratic in q and quadratic in the inter- 
action, which leads to less singular contributions to the 
free energy than the one we keep. Note the frequency re- 
strictions inherent in Eq. (2.10): sgn(a;„j) = sgn(aj„2) = 
-sgn(tj„3). Here 



'(512 i 



(^) = ^E^i(p)^2(p-fc) 



V 



(2.11) 



with aj„j a;„2 < implied, and 

^f3)2(fcl,fc2) = ;^^Gi(p)G3(p-fci)G2(p-fci-fc2) 

(2.12) 
where Gi{p) = G{p,iujni) is the single-particle Green 
function. ipi2 has a scaling form 

27rG 

fi2{k) = Nf ——ipd{Gifli-2/k) 



= (y9(fc,rJi_2) 



(2.13) 



where G is a coupling constant whose bare value is the 
inverse Fermi velocity, G — l/vp, Np is the density of 
states per spin at the Fermi level, and fli-2 — ^ni ~'-^n2- 
In d — 2, 3, and for free electrons, we find explicitly 



'Pd=2iz) = sgn(Imz)/Vl-z2 , (2.14a) 



^d=3{z) 



In 



1-2 



(2.14b) 



2 "'\-l-z^ 
which we recognize as the hydrodynamic part of the Lind- 



hard function. Equations (2.131 and (2.141 reflect the 



soft particle-hole excitations with a linear momentum- 
frequency relation in a metallic electron system. In par- 
ticular, (/3(fc, J7„ = 0) ex 1/|A;|, and ip{k = 0,rj„) c>c 
l/rin-^^ For later reference we also note the following 
identities that hold for a special form of (/p^'^-' : 



(pf^j_{k,-k) 



-'P2i2ik..-k) 



d 



diuj„ 



Vi2ik) 



(2.15) 



The fermionic action can be expressed in terms of q 
and P, and by using Eq. (2.10) and its generalizations 



to higher order one obtains a fermionic soft-mode action 
entire in terms of q. For our purposes we need only the 



Gaussian part of this action, which reads 

3 

AAq] = '8J2T. E E^9i2W ri2,34(fc);'734(-fc) ■ 

k 1.2 r=0,3 i=0 

3,4 

(2.16a) 
Here 1 = ni etc., and the Gaussian vertex is given by 



rk34(fc) = ¥'r2'(fe) + '5l-2,3-4 2T7^ 



(2.16b) 



with 7*^^^ — —7s and 7'=i'2,3 _ ^^^^^ where 7s > and 
jt,i > are the spin-singlet and spin-triplet interaction 
amplitudes. The fermionic Gaussian propagator is given 
by the inverse of the vertex. One finds 



(r9l2(fe)i934(-'«)) = Y^'^'-s'J.y 



Sl3S24Vl2ik) - 27'T5i_2,3-4 



yi2(fc)'^34(fc) 

i-2rx^^l,{k) 



where 



X^°l^ik) ^ X^"Hfe,^l-2) = -TY,Sl-2,3-4V34ik) 



(2.17a) 



(2.17b) 



34 



We see that the q-propagator is given in terms of ip, and hence is soft. The fields ^ that enter P, Eq. (2.10), are 
characterized by Gaussian propagators 



and 



{U,,ik)iq,,i-k)) - (;9i2(fc)i^34(-fc)) = Cq,2ik)iq34i-k)) 



-1 



{lA2ik)U34i-k)) = ^l'TS,^2,3 



i-2rxi°J2(fe) 



(2.17c) 



(2.17d) 



The last expression is just the interacting part of the 
g-propagator, Eq. (2.17a), as was mentioned after Eq. 
pJbl). 



Eq. (2.4b). The other possibility, which is conceptually 



The interaction amplitudes in the Gaussian fermionic 
vertex, Eq. (2. 16b I, warrant some comments. First, we 
note that the three spin-triplet amplitudes 7^. ' ' are in 
general not identical in a cyrstalline solid, and they do 
not need to be for what follows. Second, we comment 
on the two cases that result from the magnetism being 
caused by the conduction electrons, or by electrons in 
a band different from the conduction band, respectively. 
Let us first assume the latter case, which is the concep- 
tually more straightforward one. Then .Apig], which de- 
scribes the conduction electrons, is independent of the 
magnetism and contains interactions in both the spin- 
singlet and spin-triplet channels. The only restriction is 
that the latter are weak enough to not lead to magnetism 
by themselves. The conduction electrons are affected by 
the magnetization, which acts as an effective magnetic 
field, and this is described by the Zeeman coupling term. 



more complex, is that the magnetism is caused by the 
conduction electrons themselves. In this case the mag- 
netic order parameter and the soft modes q describe de- 
grees of freedom for electrons in the same band. The 
magnetic order parameter then should be thought of as 
deriving from the spin-triplet interaction between the 
conduction electrons, e.g., via a Hubbard-Stratonovich 
decoupling of the latter. This leaves the bare action Af 
with a spin-singlet interaction only. However, as long as 
the latter is present, a spin-triplet interaction will always 
be generated under renormalization. The action Af will 
therefore again contain a spin-triplet interaction ampli- 
tude, albeit one that is much weaker than the one in the 
underlying action that describes the system before the 
separation of magnetic and fermionic degrees of freedom. 
This is the case that was discussed, for ferromagnetism, 
in Ref. I13[ which used phenomenological and symmetry 
arguments to construct the fermionic part of the action. 
Finally, we mention that we assume the conduction elec- 



trons, in the absence of a nonzero magnetization (i.e., 
with the couphng constant c in Eq. (2.4b I put equal to 



zero), to indeed have three soft spin-triplet excitations at 
T = 0, which are given by Eqs. (12. 171) with i = 1,2,3. 



This is not necessarily the case. For instance, an external 
magnetic field gives two of these three channels (the ones 
transverse to the field) a mass, and a small concentration 
of magnetic impurities will make all three channels mas- 
sive without having significant other effects. However, 
in general the energy scales associated with these effects 
will be small, and they will lead to a small reduction, but 
not a complete suppression, of the tricritical temperature 
in Fig.[l] We will discuss this point in more detail in Sec. 

nil 



C. Coupling between the order parameter and the 
fermionic soft modes 



D. Generalized Mean-Field Theory 

An effective action, ^cff [-W] in terms of the order pa- 
rameter alone can be obtained by integrating out the 
fields q, 



Aes[M] = In f D[q] e-^LGw[M,9l _ 



(2.20) 



In general the evaluation of this expression is very diffi- 
cult. However, it can be evaluated exactly within a gener- 
alized mean-field approximation that was first employed 
in the context of liquid crystals and superconductorgSl 
and is defined as follows. First, we ignore temporal and 
spatial variations of the order parameter, i.e. we treat 



the fields m{x) and n{x) in Eq. (2.3) as numbers. If we 



assume ordering in the 3-direction, we have 



We are now in a position to separate the Zeeman term, 
Eq. (2.4b), into parts where the order parameter couples 



to soft and massive fermionic modes, respectively. If we 
define a temporal Fourier transform of the magnetization 
field m by 



mn{x) 



T I dT e*""^ m{x, t) 





(2.18) 



with f2„ = 2T:Tn a bosonic Matsubara frequency, then 



we can write Eq. (2.4b) in the form 



3 

Ac[M;Q] = 2cVf dxY.Y.'^nix) 
X I^(-l)''^'I]tr[(r,®s,)Qm,m+nN] -(2.19) 



r=0,3 



By expressing Q in terms of q and P by means of Eq. 
(2.8), and P in terms of ^ by means of Eq. (2.10), we 



obtain the desired coupling Ac[M, q] between the order- 
parameter fiuctuations and the fermionic soft modes q. 



MHx) 



which implies 



3j3 



N 



cos(fej • x) 



(2.21a) 



m^(a;) w Si^ (5„o ■m/Vf . 



(2.21b) 



This mean-field approximation for the order parameter 
means that only the part of Q that is diagonal in fre- 
quency space, i.e., Pmim contributes to Eq. (2.19). This 



in turn means that the contribution to P that is linear 
in g, which we had dropped from Eq. (2.10), does not 



contribute. Second, we restrict ourselves to quadratic 
order in q. That is, we treat the fermionic soft modes in 
a Gaussian approximation with a fixed magnetic order 
parameter. The validity of these approximations will be 
discussed in Sec. IIIIBI 

With these approximations the action Ac that couples 
q and the order parameter is quadratic in q and can be 
written 



Here 



AKg] = 8 Y, ^;gi2(fc);ir;2,34(fc)i934(-fe) 

r,s— 0,3 i,j 





M I 1 

0-100 




rir;2,34(fe) = ^i3'524 4cmf "^ ^ 1 I Q -1 I "^m C^. -Z^) V12' W : 



(2.22a) 



(2.22b) 



and we have used Eq. (2.15 ). The matrices give the values of J^F for the 4 possible values of (r, s) and the 16 possible 
values of (?, j). 



The integral over q in Eq. (2.20) can now easily be carried out. For the free-energy density / = —TAcs/y we 
obtain 



/ = fo{m,n) + Af{m) 



(2.23a) 



Here /o = —TAop/V is the mean-field free energy in the absence of a coupling to the fermionic soft modes. For 
Af{m), which is the contribution to the free energy due to this coupling, one finds 



A/(m) = :^^V^lniV(fc,f^„;m) 



(2.23b) 



where J^k denotes a wave vector sum such that |fe| < A with A an ultraviolet cutoff, and 

Nik,n^;m)^-l6c''^tA-ft,2m^nl[ip<-''\k,n„))\-Hk,^n) + ^-^{k,n„)ll [l-27t,,x^°)(fc,f^„)] . (2.23c) 



The equation of state is obtained by minimizing the 
free energy density. In the absence of a coupling be- 
tween the order parameter and the fermionic soft modes 
this amounts to minimizing /g , which yields the ordinary 
mean-field equation of state. For a ferromagnet, the lat- 
ter has the usual Landau form. For a ferrimagnet, the 
equation of state depends on details of the magnetic or- 
der. It can be complicated and describe several different 
phases, see, e.g., Ref. J^ However, generically the first 
phase encountered as one approaches from the paramag- 
netic state is entered via a second-order transition. After 
minimizing /o and expressing n in terms of m one thus 
has again an ordinary mean-field equation of state given 
by 



0{m^) 



(2.24) 



where h is an external magnetic field in the 3-direction, 
w > 0, and the transition occurs at r = O.l^SIln Appendix 
[X]we recall a very simple model that leads to this result. 
The second term on the right-hand side of Eq. (2.23a) 



gives an additional contribution to the equation of state, 
which then reads 



h — r m + um — 64:171 c jt,i'~ft,2 



V 



k 



N{k,nn;m) 



(2.25) 



This is the desired generalized mean-field equation of 
state which takes into account the coupling of the order 
parameter to the fermionic soft modes. 



E. 



Discussion of the Generalized Mean-Field 
Equation of State 



With some effort the integrals in Eqs. (2.23b I and 



(2.25) can be explicitly performed. However, the salient 



points can be seen by simple scaling considerations and 
dimensional analysis. Equations (2.11) and (2.13) im- 



ply that the frequency 0„ scales as the wavenumbcr 
/c, J7„ ~ k, and that </?(/?, $!„) ^ 1/fc ^ 1/f^n, which 
also can be seen explicitly from Eqs. (|2.14|. Equation 



( |2.15) impl ies that (p^^'^ik,^^) - 1/P - l/0^_ Equa- 
tion ( |2.23c ) then shows that there is a length scale L™, 
or a corresponding frequency scale ujm, that scales as 
Lm ^ l/i^m ^ 1/m. If one attempts to expand A/(m), 
Eq. (2.23b), in powers of m at T = 0, then nonanalytici- 



ties will occur at next-to-leading order for all d < 3. 

An alternative way to describe this mechanism is to 
say that of the three soft fermionic spin-triplet excita- 
tions, Eq. (2.17a) with r = s = 0, 3 and i = j = 1,2,3, 



two (namely, the ones transverse to the order parameter 
direction) acquire a mass due to the coupling between 
the fermions and the order parameter to, as can be seen 
explicitly from Eq. (2.22b). This acquisition of a mass 



by a generic soft mode due the spontaneous breaking of 
a continuous symmetry is an example of the Anderson- 
Higgs mechanism ]^ "^ even though the broken symmetry 
in this case is not a gauge symmetry, see the discussion 
in Sec. |III A| It implies in turn that the free energy is a 
nonanalytic function of to,. 

At nonzero temperatures the singularities are cut off 
by T according to to. ^ T. That is, a crossover occurs 
from TTi-scaling to T-scaling when the Zeeman splitting 
is comparable to the temperature, or the thermal length 
scale Lt oc 1/T is comparable to the magnetic length 
scale Lfn mentioned above. Taking into account the sign 
of TV, Eq. (2.23c), one finds schematically, for 1 < d < 3, 



A/(TO) = -VTO2(TO2-Fr2)('^-l)/2 , 

and for d = 3 

Af{m) = — TO \n{m + T") , 



(2.26a) 
(2.26b) 



where w > is a positive constant. 

The most important aspects of this result, as far as 
the order of the transition is concerned, are the sign of 
V and the power of m at T = 0. For all d < 3 there 
is a negative term in the free energy that dominates the 
m^ in the Landau free energy and hence necessarily leads 
to a first-order transition. Another way to see this is by 
expanding A/(to,), Eq. (2.26a), in powers of tti for T > 0. 



The leading term is proportional to —w^/T^~'^. That is, 
there is a negative to term whose prefactor diverges as 
T — )■ for all d < 3, which implies that there will be a 
tricritical point at some temperature. The free energy for 








m 



FIG. 2: Schematic sketch of the free energy for three values 
of the parameter r. The first-order transition occurs at r = 
ri > 0. It pre-empts the second-order transition of Landau 
theory which would occur at r = 0. 



magnetic impurities. We note in passing that an inter- 
esting system is provided by the easy-plane ferromagnet 
URhGe, where an in-plane magnetic field transverse to 
the magnetization has been used to tune the transition, 
access the tricritical point, and map out the tricritical 
wings.'^ This situation requires a refinement of the the- 
ory presented above, which will be reported elsewhere .1211 



III. DISCUSSION, AND CONCLUSION 

We now discuss our results, before concluding with a 
summary. 



three different values of r is plotted schematically in Fig. 
[2] For this schematic free energy, the equation of state 
in the case d = 3, for which many experimental results 
exist, takes the form 



The mechanism behind the first-order transition 



h = rm-\ — m ln(m -|- T ) 



m M H — 



4 m2 + T2 



{d = 3) (2.27) 



Also of interest is the other physical dimensionality, d 
2, where the equation of state reads 



rm — 2vm{m +T ) ' 



+m'' u - 



(^2+ 2^2)1/2 



(d = 2) (2.28) 



Here the analyticity is stronger than in the ?>-d case, with 
a negative rn^-term in the equation of state at T = 0. 
This is particularly interesting in the case of Ising mag- 
nets, which display long-range order in d = 2 even at 
T > 0. The case of Heisenberg and XY magnets, which 
do not show true long-range order in d = 2 except at 
T = 0, is more complicated. 

These are the same results that were obtained using 
a more phenomenological theory of the fermionic soft 
modes in Ref.[TSl They were discussed extensively in that 
reference, as well as in Refs.|4]and[7l There is no need to 
repeat this discussion here, and the salient features are 
summarized by the schematic phase diagram shown in 
Fig. [Tl The important conclusion of the current paper is 
that the validity of these results, in addition to itinerant 
Heisenberg ferromagnets, extends to metallic ferromag- 
nets where the magnetism is not due to the conduction 
electrons, to metallic ferromagnets in the XY or Ising 
universality class, and also to metallic ferrimagncts. The 
only condition is that the conduction electrons are not 
subject to strong spin-symmetry breaking effects such as 



leads to the first-order tran- 
II E is precisely analogous to 



The mechanism that 
sition discussed in Sec. 
the fluctuation-induced first-order transition discussed in 
Ref. for the BCS-superconductor transition and the 
nematic-to-smectic-A transition in liquid crystals. An 
important physical ingredient is an underlying "generic" 
soft mode, i.e., one that is not related to the phase tran- 
sition in question, but couples to the order parameter. 
In the case of liquid crystals this soft mode is the ne- 
matic Goldstone mode, in the case of superconductors, 
the vector potential, in the present case, the spin-triplet 
particle-hole excitation. At the transition of interest, this 
soft mode acquires a mass that is given in terms of the 
nonzero expectation value of the order parameter. This 
general mass-generating mechanism was first pointed out 
by Anderson, and is now known as the Anderson-Higgs 
mechanism.'SSHlll This coupling of the order parameter to 
underlying soft modes leads to a non-analytic term in the 
Landau free energy that is dominant over the usual quar- 
tic term and has a negative sign, leading to a first-order 
transition. It should be stressed that this is only one way 
to realize a fluctuation-induced first-order transition; an- 
other one, for instance, is realized by a 0^-theory with a 
cubic anisotropyl^ The current realization is analogous 
to the case of scalar electrodynamics studied by Cole- 
man and Weinberg in a particle-physics context .'2S1 It is 
also worthwhile noting that the analogy between super- 
conductors on one hand, and liquid crystals and quantum 
magnets on the other, breaks down in the ordered phase. 
In the former case, the Goldstone mode gets absorbed 
into the longitudinal component of the vector potential, 
which is massive, and there is no soft mode in the ordered 
phase. In the latter, there are Goldstone modes in the 
ordered phases, namely, a "smecton" with an anisotropic 
dispersion relation in the smectic-A phase (Ref. |271 see 
also Ref. .28() and magnons in the magnetic phase. 
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B. Universality of the first-order transition, and 
the validity of the generalized mean-fleld theory 

Experimentally, all examples of clean low-Tc ferromag- 
nets (for disordered systems, see below; ferrimagnets so 
far have not been systematically studied from this point 
of view) show a first-order transition if the Curie temper- 
ature is suppressed far enough. There is not a single ex- 
ample of a quantum critical point in zero magnetic field. 
While this is consistent with the generalized mean-field 
theory theory presented in Sec. ITT] it is somewhat surpris- 
ing when compared with the case of liquid crystals, where 
an analogous theory also predicts a first-order transition. 
In this case, in stark contrast to that of quantum mag- 
nets, the observed transition is usually of second order, 
and only recently have examples of a (weakly) first-order 
transition been found.'^Sl These observations beg the ques- 
tion whether in the case of quantum magnets the gener- 
alized mean-field approximation is more generally valid 
than in classical systems. 

To discuss this point, we first observe that we have 
made three approximations to treat the action given by 
Eq. (2.1a|. First, we have integrated out the fermionic 



massive modes in a saddle-point approximation that re- 
spects the Ward identi ty th at governs the soft-mode 
structure of the system.'^S'^ Second, we have kept the 
soft fermionic degrees of freedom only to Gaussian or- 
der in the soft modes q. Third, we have treated the 
order parameter in a mean-field approximation. These 
approximations are not independent of one another, and 
the first two simplifications do not constitute any addi- 
tional approximation over and above the last one. This 
can be seen as follows. 

The mean-fleld approximation for the order parame- 
ter means that the fermionic degrees of freedom describe 
an interacting electron system that is spin-polarized by 
the coupling to the homogeneous magnetization, which 
acts as an effective external magnetic field. The state 
of the fermionic subsystem is thus described by a stable 
Fermi-liquid fixed point. Corrections to the fermionic 
soft-mode action due to massive degrees of freedom are 
irrelevant with respect to this fixed point by at least one- 
half power of frequency or wavenumber in all dimensions, 
and thus cannot change the properties of system.'^^ Sim- 
ilarly, only the terms quadratic in q contribute to the 
fixed-point action; all higher-order terms are irrelevant 
by power counting. Keeping terms of higher order in q 
will therefore renormalize the parameters of the theory, 
but it cannot change its structure. In particular, it can- 
not change the sign of the term in the equation of state, 
Eqs. ( 2.27[ 2.28), that is due to the soft fermionic fluctu- 
ations and leads to the first-order transition. 

This leaves the mean-field approximation for the order 
parameter to be discussed. If the first-order transition 
at r = ri occurs far from the second-order transition at 
r = that is pre-empted by it (see Fig. [2]) , then order- 
parameter fluctuations are negligible and the results of 
the generalized mean-field theory are qualitatively cor- 



rect. If, however, the first-order transition occurs close 
to the putative second-order one, i.e., if the minimum 
in the free energy in Fig. [2] is very shallow, then it is 
less clear whether order-parameter fluctuations can be 
neglected. One key difference between classical liquid 
crystals and quantum magnets is that in the former case, 
the system is below the upper critical dimension d^ = 4 
for the (unrealized) phase transition that would occur 
in the absence of any coupling between the smectic or- 
der parameter and the nematic soft modes. In contrast, 
the quantum magnetic systems are above the correspond- 
ing upper critical dimension d^ = 1 that follows from 
Hertz theory, and even with that coupling taken into ac- 
count, ordinary mean-fleld theory becomes exact, as far 
as the description of the phase transition is concerned, 
for d > S.Ei This strongly suggests that order-parameter 
fluctuations are of much less importance in the case of 
quantum magnets, and it provides a possible explanation 
of the fact that the observed transition is universally of 
flrst order. 

Irrespective of these observations, the role of order- 
parameter fluctuations in quantum magnets is a topic 
that warrants additional work. For the case where the 
magnetism is not produced by the conduction electrons, 
this will require an action that properly describes lo- 
calized magnetic moments and their fluctuations, e.g., 
the one given in Ref. 1321 For itinerant magnets, i.e., if 
the magnetism is due to the conduction electrons them- 
selves, the theory developed in Sec. [IT] will apply, but 
the order-parameter fluctuations and the fermionic ex- 
citations both need to be kept, along the lines of the 
phenomenological theory of Ref. [131 The latter reference 
gave a scenario that can lead to a second-order transi- 
tion in the magnetic case. It would also be interesting to 
experimentally study quantum ferromagnets or ferrimag- 
nets in d = 2, where order-parameter fluctuations will be 
stronger than in rf = 3. 



C. The effects of quenched disorder 

So far we have discussed the case of clean or pure mag- 
nets. Impurities, modeled by quenched disorder, have 
important effects that are both needed to understand ex- 
perimental observations in certain systems, and to pre- 
dict effects that can serve to ascertain that the flrst-order 
transition in pure samples is indeed due to the posited 
mechanism. 

Quenched disorder changes the soft-mode spectrum of 
the fermions. It gives the ballistic soft modes that are 
represented by Eqs. (2.17) as mass, and leads to new soft 
modes that are diffusive. In the context of the current 
theory, this change has two principal effects. First, it 



cuts off the nonanalyticity in the clean equation of state, 
Eqs. (2.27 2.28). Second, it leads to a new nonanalytic 



term in the equation of state that has the opposite sign 
and whose prefactor vanishes in the clean limit .l^iJ The 
resulting schematic generalized Landau theory has been 
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discussed in Ref. |31 A more detailed model discussion 
that allows for semi-quantitative predictions of the effects 
of disorder will be presented elsewherepil here we just 
present the most pertinent aspects of such a model cal- 
culation. A good representation of the mean-field equa- 
tion of state for realistic values of the magnetization, the 
temperature, and the disorder, is 

?! / m.^ 



+ -m^ \n[cm'^ + {l/kFe + bTf] +um^ , (3.1) 



which generalizes Eq. (2.27) in the presence of quenched 



disorder. Here the magnetic field h and the temperature 
T are measured in units of the Fermi energy ep and the 
Fermi temperature Tp, respectively, and the magnetiza- 
tion m is measured in units of the conduction electron 
density (we put /ie = 1). The dimcnsionlcss coupling 
constant v is proportional to the fourth power of the ef- 
fective spin-triplet interaction amplitude of the conduc- 
tion electrons. It is a measure of how strongly corre- 
lated the conduction electrons are, and it is bounded 
above by a stability criterion that requires v < 0.5. 
k-p is the Fermi wave number of the conduction elec- 
trons, and £ is the elastic mean-free path. Within a 
Drude model, and for good metals, one has approxi- 
mately kpE « 1,000/ {pQ/^ricm), with po the residual 
electrical resistivity, c and b are dimensionless constants 
that are equal to c = 1/45 and 6 = Stt in a model 
calculation.^^ The second factor in the second term on 
the right-hand side is a reasonable representation, for re- 
alistic parameter values, of a more complicated scaling 
function 



m ' g{kp£m,bT/m) 



1,3/2 



(fer)3/2 



(3.2) 



that depends on the disorder in addition to the temper- 
ature, and wc have dropped the last term in Eq. (2.27) 
from Eq. (3.1 ) since one generically expects v <^ u. 
At T = 0, and in a clean system, Eq. (3.1 1 yields a 



first-order transition at ri = wto^/4, where the magneti- 
zation discontinuously jumps from m ^ to m ^ nii = 
g_(i+2„/„)/2 -^j^j^ y _ QJ4 ^^^ ^ _ Q Q2 ^-his yields 

nil ~ 4 X 10^'^, which is reasonable for a weak ferromag- 
net. Similarly, there is a tricritical temperature given 
by Ttc/Tp = (1/6) exp(— u/w); with the same parame- 
ter values this yields Ttc/Tp « 10"^, or Ttc « lOK for 
Tp = 100, 000 K, which is also reasonable. This tricritical 
point gets destroyed by quenched disorder on the order of 
fcp£ w hTtc/T-p K, 1,000, or a residual resisitivity on the 
order of po w 1/iilcm. At this point the second term on 
the right-hand side of Eq. ( 3.1 ) is still very small, and the 



critical behavior at the resulting quantum critical point 
is given by ordinary mean-field exponents except extrely 
close to the transition, where it crosses over to the crit- 
ical behavior derived in Ref. [TT] For instance, in this 
asymptotic region the critical exponents /3 and 5, defined 



by m{h = 0) ex \r\^ and m{r = 0) ex h}/^ , respectively, 
are given by /3 = 1/2 and 5 = 3/2, as opposed to the 
mean-field values /3 = 1/2 and (5 = 3. Only for sub- 
stantially larger values of the disorder, po ~ lOOpilcm 
with the above parameters, does the asymptotic critical 
behavior extend over a sizeable range of r values (up to 
\r\ w 0.01). This observation explains why an experi- 
ment on Nia;Pdi_a:, which shows a ferromagnetic transi- 
tion at a very small value oi x {x ~ 0.025) corresponding 
to weak disorder, found mean-field exponents consistent 
with Hertz theory,E3 whereas Bauer et al.SSl found non- 
mean-ficld exponents, at least some of which were con- 
sistent with Ref. [T71 in URu2-£!;Rea;Si2, where the ferro- 
magnetic transition occurs at x « 0.15 with the residual 
resistivity on the order of po ~ lOO/iricm.lSSl 



D. Conclusion 

In conclusion, we have extended a previous theory 
of quantum ferromagnets in several important ways. 
We have shown that the mechanism that leads to the 
paramagnet-to-ferromagnet transition at low tempera- 
ture in d = 3 and d = 2 to be generically of first order, 
which was first reported in Ref. HI is valid in anisotropic 
ferromagnets, in ferrimagnets, and in metallic ferromag- 
nets where the conduction electrons are not the source of 
the magnetization, in addition to the case of isotropic 
itinerant ferromagnets originally considered. This ex- 
plains why the low-temperature transition is observed 
to be of first order in highly anisotropic ferromagnets, 
and it much expands the class of materials for which this 
phenomenon is predicted. For clean magnets, an effec- 
tive theory of soft fcrmionic modes recently developed 
in Ref. il6i, has provided a technical basis that improves 
on the phenomenological theory of Ref. [131 In the pres- 
ence of quenched disorder, the theory allows for a semi- 
quantitative description of the suppression and ultimate 
destruction of the tricritical point. A sizeable range of 
disorder exists where the observable critical behavior is 
predicted to be mean-field like, whereas for very large 
disorder the asymptotic critical region, which is char- 
acterized by non-mean-field Gaussian critical exponents, 
expands and eventually eliminates the mean-field region. 



Appendix A: A simple mean-fleld model of a 
ferrimagnet 

Here we recall a very simple mean-field model of the 
transition from a paramagnet to long-range ferrimagnetic 
order.lJ^ Consider a one-dimensional chain of alternating 
magnetic moments /x^, /i^ that are antiferromagnetically 
coupled. Weiss theory assumes that the a-moments and 
^-moments are subject to effective magnetic fields 



Ba 

Bb 



-XMb 

'XMa 



(Ala) 
(Alb) 
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respectively, where A > 0. The magnetizations Ma^b are 
given by the BriUouin expressions 



Ma = 1^ fJ-a tanh(^a H/T 
Mh = V fJ-b tanh(/ih H/T - 



- fia Ba/T) , (A2a) 
fib Bb/T) . (A2b) 



Here H is an external magnetic field, T is the tem- 
perature, and V is the number of magnetic moments 
of each species. If one defines reduced magnetic fields 
ha,b = H/v^Jia.b^, & reduccd temperature t = T /v^a^^b^, 
and reduced moments nia^b — Ma,bli^J£aj}, then one sees 



that the Weiss mean-field equations (Al A2) have a so- 
lution nia = —mb = in, where m is the solution of the 
usual mean-field equation of state 



rm + fh'^/i + 0{rrf) 



(A3) 



where r = t — 1. This simple model thus describes a 
transition at t = 1 to ferrimagnetic order where the 
homogeneous magnetization is given by m = Ma + 
Mb = h'ifJ-a — tJ'b)'iTi and the staggered magnetization 
n = Ma — Mb — v(pa + Mh)"^i is proportional to m. 
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